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1. Introduction 

One of the most promising candidates for quantum computation is the implementation 
using spin systems as quantum bits (qubits) [H |2l [3l IH O |6l [7]. Combining with 
nanostructure technology, they have the potential advantage of being scalable to a 
large system. In particular, the spin in a quantum dot exhibits a relatively long 
coherence time compared with fast gate operation times. This makes it a good 
candidate for a quantum information carrier [HI [9l HOl [H]. However the influence of 
the environment, especially the spin environment, on a spin, which usually causes the 
decoherence of the spin, is inevitable. Typically, the material environment must be 
present in order to host the spin or to locally control the electric or magnetic fields 
experienced by the spin. Therefore the decoherence behavior of a central spin or 
several spins interacting with a spin bath has attracted much attention in recent years 

The problem of a central spin coupled to an antiferromagnetic (AF) spin 
environment was investigated in Refs. [261 EH I2H]- In these investigations, a pure 
dephasing model in which the self-Hamiltonian of the undriven central spin commutes 
with the interaction Hamiltonian to the environment was considered. A similar problem 
of a spin- 1 impurity embedded in an AF environment was studied in the context of 
quantum frustration of decoherence in Ref. [21]. There the impurity spin is coupled 
locally in real space to just one spin of the AF environment, in contrast to the central 
spin model [3nil3ll|3a[33l[3a[35l|36l|371[Ml[39ll^ where the central spin is coupled 
isotropically with equal strength to all the spins of the environment. In this paper, we 
consider a more general case of a central spin (qubit) driven by an external field and 
coupled to an AF spin bath (environment). To control the quantum states of a qubit, 
some types of time- dependent controllable manipulations , which can be electrical [42] , 
optical [13], or magnetic [B], are desirable. Thus it is important and necessary to 
study the behavior of a qubit (central spin) or several qubits (spins) interacting with an 
environment (a spin bath) in the presence of a driving field [45} |46l |47l HH] . Recently, the 
coherent dynamics of a single central spin (a nitrogen-vacancy center) coupled to a bath 
of spins (nitrogen impurities) in diamond was studied experimentally [47J. To realize the 
manipulation of the central spin, the pulsed radio-frequency radiation was used. Under 
the control of a time-dependent magnetic field, Ref. [H] considered a central spin system 
coupled to a spin bath. The decay of Rabi oscillations and the loss of entanglement were 
discussed. However, the correlations between spins in the spin bath were neglected. 
Recent investigations indicated that the internal dynamics of the spin bath could be 
crucial to the decoherence of the central spin [ia[ia[Ta[T5l[T6l[I71ll8l[l9l|2nil2lll22]. In 
our driven spin model, the interactions between constituent spins of the AF environment 
are taken into account. 

The self-Hamiltonian of the central spin, after a transformation to a frame rotating 
with the frequency of the driving field, can be written in a form of Hs = eSq+qSq in the 
rotating wave approximation. Due to the additional driving term of qSq that provides 
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energy into the system and does not commute with the interaction Hamiltonian that 
couples to the AF environment in our model, the dynamics of the central spin in this 
case is dramatically different from the undriven pure dephasing behaviors investigated 
in Refs. [261 HZ] • It was shown in Ref . [17] that the form and the rate of Rabi oscillation 
decay are useful in experimentally determining the intrabath coupling strength for a 
broad class of solid-state systems [T7]. This is also the case in our problem. After the 
use of the spin- wave approximation to deal with the AF environment, we employed an 
elegant mathematic technique to trace over the AF environmental degrees of freedom 
exactly to obtain the reduced density matrix of the driven spin. This enables us to 
study and describe the decay behaviors of the Rabi oscillations for different initial 
states and parameters of the central spin-AF environment model beyond the Markovian 
approximation and Born approximation (perturbation). With the reduced density 
matrix obtained nonperturbatively, we also investigate and discuss the decoherence and 
the pointer states of the central spin from the perspective of the von Neumann entropy. 
We find that the decoherence behavior of the central spin depends on the detuning, 
driving strength, the coupling between the central spin and the spin environment, and 
an important factor associated with the number of the spins, the lattice structure 
and the temperature of the environment. If the detuning exists, the Rabi oscillations 
may show a behavior of collapses and revivals; however, if the detuning is zero, such 
a behavior will not appear. We investigate the weighted frequency distribution of the 
time evolution of the central spin inversion and give this phenomenon a reasonable 
explanation. The form and the rate of Rabi oscillation decay is useful in determining 
the intrabath coupling strength and other related properties of the qubit-environment 
system [TTj. Although we concentrate on the central spin model, our study is applicable 
to similar models of a pseudo spin or a qubit coupled to an environment. For example, in 
the study of the decoherence behavior of a flux qubit interacting with a spin environment, 
the same self-Hamiltonian of the qubit, Hs = eSq+qSq, can be used [25]. In such case, e 
is the bias energy and g is tunnel splitting of the flux qubit, and the two eigenstates of 
correspond to macroscopically distinct states that have a clockwise or an anticlockwise 
circulating current [25] which can be denoted as |1) or |0). To make contact of the 
driven central spin model with experiments, one may envisage a setup of a small ring- 
shaped flux qubit located at a distance above or below an also ring-shaped AF material 
that has a common symmetric z-axis with the flux qubit. In this case, the coupling 
strength between the flux qubit and each of the constituent spins of the AF material 
(environment) may be regarded to be almost the same. 

The paper is organized as follows. In Sec. [21 the model Hamiltonian is introduced 
and the spin wave approximation is applied to map the spin operators of the AF 
environment to bosonic operators. After tracing over the environmental modes, the 
reduced density matrix is obtained and the dynamics of the central spin is calculated. 
We also investigate the decoherence and the pointer states of the central spin for the 
cases of zero detuning and nonzero detuning, and calculate the von Neumann entropy 
which is a measure of the purity of the mixed state. Numerical results and discussions 
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are presented in Sec. [21 Conclusions are given in Sec. ID 



2. Model and Calculations 

2.1. Model and transformed Hamiltonian 

We consider a central spin driven by an external microwave magnetic field and embedded 
in an AF material. To detect the central spin and control its states, the frequency Uc of 
the microwave magnetic field is tuned to be resonant or near resonant with the central 
spin. Furthermore, the central spin and the AF environment are made of spin-i atoms. 
The total Hamiltonian of our model can be written as 

H = Hs + HsB + Hb, (1) 

where Hs, Hb are the Hamiltonians of the central spin and the AF environment 
respectively, and Hsb is the interaction between them [26l l27l |49| [50] . They can be 
written as {h= 1) 

Hs =/io5o+^7(5o^e--=* + 5o-e*-=*), (2) 
ji 

Hsb = - 77^*^0 X^('S'a,i + Sb,i)y (3) 

Hb = Sa,i ■ + ■ ^a,j+5 ' (4) 

where /xq is the Larmor frequency and represents the coupling constant with a local 
magnetic field in the z direction. The magnetic field creates a local Zeemann splitting, 
which can then be accessed by the driving field with frequency Uc and (real) coupling 
strength g. The coupling strength g is proportional to the amplitude of the driving 
field. We assume that the spin structure of the AF environment may be divided into 
two interpenetrating sublattices a and h with the property that all nearest neighbors 
of an atom on a lie on 6, and vice versa [51j. Sa,i {^b,j) represents the spin operator 
of the iih (jth) atom on sublattice a (b). The indices i and j label the N atoms in 
each sublattice, whereas the vectors 5 connect atom i or j with its nearest neighbors. 
J is the exchange interaction and is positive for AF environment. The effects of the 
next nearest-neighbor interactions are neglected. For simplicity, significant interaction 
between the central spin and the environment is assumed to be of the Ising type. This 
type of interaction has gained additional importance because of its relevance to quantum 
information processing [211 [l6l HHl [52]. The coupling constant between the central spin 
and AF environment is scaled as Jq/\/N such that nontrivial finite limit of — t- oo 
can exist [26l [271 SHI ESI [53] . For convenience, we denote in the following the scaled 
interaction between the central spin and the AF environment as Jq = Jq/^/N. 

In a frame rotating with the frequency of the driving field, the Hamiltonian of the 
central spin becomes 

Hs = eS^ + gS^,, (5) 
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where the detuning 
e = /io 

Using the Holstein-Primakoff transformation 

St. = 



1 



s: 



St, = ^^i-btb„ S~^ = ^1 - hlb,b,, 



S'„ 



--aiai 



(6) 
(7) 

~ "j \l ^ "7 "J' "^6,7 ~ \/ ":i "J "J' '~^b,j ~ "7 "3 o' 

we map spin operators of the AF environment onto bosonic operators. We will consider 
the situation that the environment is in the low-temperature and low-excitation limit 
such that the spin operators in Eqs. ([7]) and ([8]) can be approximated as S^^ ~ Oj, 
and S^j ~ . This can be justified because in this limit, the number of excitation is 
small, and the thermal averages < a^Oj > and < blhi > are expected to be of the order 
0{1/N) and can be safely neglected when is very large. The Hamiltonians Hsb and 
Hb can then be written in the spin-wave approximation [51] as 



i 

Hb = ~ \nMJ + MJY,{(^t(^^ + Kh) + J Y.{a,h^+i + a+b 



(9) 

(10) 



i,6 



where M is the number of the nearest neighbors of an atom. We note here that in 
obtaining Hamiltonian f lTOj) in line with the approximations of 5*^^ ~ Oj, and S^j ~ in 
the low excitation limit, we have neglected terms that contain products of four operators. 
The low excitations correspond to low temperatures, T <^ T/v, where T/v is the Neel 
temperature [5^. Then transforming Eqs. Qand ffTOl) to the momentum space, we have 

(11) 

k 



Hsb = - JqSq ^(&k&k - ctkCtk), 

k 

Hb = - ^NMJ + MjJ2K(^i^ + Kbi. 

k 

+ MJ ^ 7k (a+ 6+ + ak6k), 



where 7k 



where 



M ^ '^s^'^^^- Furthermore, by using the Bogoliubov transformation, 

«k = ^ikCtk - fk&k' 

/3k = Mk&k - ^^k^k , 
(1 + A)/2, vl = 



(12) 



(13) 
(14) 



-(1 — A)/2, A = 1/a/1 — 7k, and the Hamiltonians Hsb 
and Hb can be diagonalized {h = 1) and be written as 



Hsb = - JoS'o ' "k"k), 

k 

Hb = - \nMJ + 5^ a;k («,^«k + /3,[/3k + l) 



(15) 
(16) 
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where aj^ (ok) and (/3k) are the creation (annihilation) operators of the two different 
magnons with wavevector k and frequency cuk respectively. For a cubic crystal system 
in the small k approximation, 

Wk = {2Mf'^Jkl, (17) 

where / is the side length of cubic primitive cell of the sublattice. The constants —^NMJ 
and can be neglected. Finally, the transformed Hamiltonians become 

Hs =eS', + gSl (18) 
HsB = - JoS^ $](/3,[/3k - «,^ak), (19) 

k 

Hb = 1] Wk («,^ak + /3,[/3k) ■ (20) 

k 

Similar to the famous spin-boson model |S3 EH], the interaction Hamiltonian does 
not commute with the self-Hamiltonian of the spin. However, a significant difference 
from the spin-boson model [551 ES] is that the interaction Hamiltonian commutes with 
the bath Hamiltonian, i.e., [Hsb,Hb] = 0. So the problem of the total transformed 
Hamiltonian, Eqs. ( fT8l) -( l20|) . can be solved exactly, even in the case of multi-environment 
modes and finite environment temperatures. 

2.2. Reduced density matrix 

We assume the initial total density matrix of the composed system is separable, i.e., 
p(0) = ps'(O) (8> Pb- The density matrix of the AF bath satisfies the Boltzmann 
distribution pB = e~^^^'^ /Z , where Z is the partition function and the Boltzmann 
constant has been set to one. If the initial state of the qubit (central spin) is taken as 
P5(0) = |V^(0)) (^(0)1 where 

|^(0)) = 5|1)+7|0), (21) 

|5p + |7r = l, (22) 

then the reduced density matrix operator of the qubit can be written as 

ps{t) = ^tiB [e-'^*|^(0))e-^-/^(V^(0)|e^^*] 

z 

+ -S-f*tTB [e-^iIis+HsB)t\l^^-HB/T^Q^^i{Hs+HsB)t^ 

z 

+ -\j\hiB [e-^(^^+^^«)*|0)e-^^/^(0|e^(^^+^«^)*] , (23) 
Z 

where tr^ denotes the partial trace taken over the Hilbert space Hb of the environment. 
The partition function Z can be evaluated as 

Z = tr^e-^-/^ 
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nr 



= e ■' ^ ) , 

where V is the volume of the environment. At a low temperature such that Ur, 



(24) 
» T, 



we may extend the upper limit of the integration to infinity. With x = {2My^'^Jkl/T 
and N = V/l^, we obtain 

^ ^ g-2n/J»ln(l-e-)x2dx^ ^^25) 



where 



n = — . (26) 

4v^7r2M3/2j3 

To obtain the exact density matrix operator, Eq. f l23|) . we need to evaluate 
Q-i{Hs+HsB)t\^\^ and e~*^^'5+^^^-'*|0). To proceed, we adopt a special operator technique 
presented in Ref. [35]. We can see that the Hamiltonian Hs + Hsb contains operators 
ttk 5 "^k, /^k , /3k5 , Sq, and Sq, where Sq and Sq change the system state from |1) to 
|0), and vice versa. It is then obvious that we can write 

^-i(Hs+HsB)t^l^ = A\l) + B\0), (27) 

where A and B are functions of operators , a^, (3^ , /3k, and time t. Using the 
Schrodinger equation identity 
.d 



dt 

and Eq. we obtain 



4^ = -B-- 

dt 2 2 



d ^ a , 1 
i—B =-^A + - 
dt 2 2 



-^o2^(/3k/3k - akttk) 



ttkttk) 



A, 



B, 



(28) 

(29) 
(30) 



with the initial conditions from Eq. (12 7p given by 

A{0) = 1, (31) 

5(0) = 0. (32) 

We note that the coefficients of Eqs. fl29|) and fl30|) involve only the operators a^a^ and 
/3,^/3k- As a result, we know that A and 5 are functions of a^a-^, P^Pk, and t. They 
therefore commute with each other. Consequently, we can treat Eqs. (l29ll and (130|1 as 
coupled complex- number differential equations and solve them in a usual way. This 
operator approach allows us to solve Eq. (1271) and obtain 



A 
B 



cos(/€t/2) + i 



Jo Ek(/5k /5k - ak^k) - £ 



K 



sin(/€t/2). 



— i— sin(Kt/2), 



(33) 
(34) 
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where 



K 



\ 



1 2 



^0 X^(/3k /3k - ak^k) 



+ 9' 



(35) 



Following the similar calculations above, we can evaluate the time evolution for the 
initial spin state of |0). Let 

^-iiHs+HsB)t\Q^ =C\l) + D|0). (36) 
In a similar way, we obtain 

C = B, (37) 
D = A+. (38) 

With Eqs. ([27D, ([33]), ([34]), ([36]), ([37]) and ([38]), the reduced density matrix operator, 
Eq. (l23l) . can be obtained analytically using a particular mathematical method to deal 
with the trace over the degrees of freedom of the thermal AF environment. This 
particular mathematical method will be described in subsection 12.31 

Note that our approach also applies to the case in which the qubit (central 
spin) is initially in a mixed state. For example, if the initial state for the qubit is 
Ps(0) = + |7p|0)(0|, the corresponding reduced density matrix is Eq. (!23ll 

provided that the second and third terms that contain, respectively, ^7* and 5*7 on its 
right-hand side are removed. 



2. 3. Dynamics of the central spin 

With the time evolution of the density matrix operator, we can investigate the dynamical 
behavior of the central spin, which is of particular interest from the perspective of 
practical application. Here we discuss the time dependence of the expectation value of 
S'o(t). It can be written as 

= (l|5oXt)|l) + (0|5oXt)|0) 

= —\6\hiB [{A+A - B+B) e-^^/^] 
2^ 

+ ^\l?^^B[{C^C-D^D)e-^-l^] 
+ ^h*tTB[{AC^-BD^) e-^-l^\ 

+ ^5*7trB {{A-^C - B+D) g-^-/^] . (39) 

As demonstrated in Eq. (12^ . by decomposing the Hamiltonian Hb and tracing over 
the modes of and (3^ separately, the partition function Z can be calculated. But 
according to the expression of Eqs. f l33]) - f[35]) . it is impossible to do so in Eq. f l39]l . 
Here we introduce a particular mathematical method to perform the trace over the 
environment degrees of freedom. Our procedure takes two steps. First, the states with 
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definite eigenvalues, say Pi and P2, of tlie respective operators Xlk'^k'^k and ^k/^k/^k 
are traced over. Tlien we sum over all possible eigenvalues Pi and P2 of the operators 
Xlk'^k'^k and Xlk/^k/^k- In this way, we obtain the final expression of (S'o(t)), i.e., 

(Sm) = ^ E E ^2) + f2{Pi. P2)] A(Pi)A(P2), (40) 



2Z 

Pl=0P2=0 



where 



/i(Pi, P2) = m'- \i\'){A*A - B*B), (41) 
/2(Pi,P2) = 4Re(57*^P*). (42) 

Here the operator «k (^k /^k /^k) in A and B has been replaced by integer Pi (P2) 

of its corresponding eigenvalues. The two conditional partition functions are defined as 

A(Pi) = trB(P,)e-S'<-''"^''/^, (43) 

A{P,) = tTBiP,)e-^-''-^-^-^^, (44) 

where only the states with eigenvalue Pi (P2) of the operator a^a]^ iJ^k /^k /^k) are 
traced over. Because of the restriction in the trace, the evaluation of A(P) is a little bit 
involved. To proceed, we define a generating function G{X) for A(P) in the following 
manner [57]. For any real number A (|A| < 1), we define 

00 

G(A) = E A''trB(p)e~ Ek-ka+a,/T_ ^45^, 

P=0 

The Taylor polynomial expansion of function G{X) at A = can be written as 

G(A) = f:A-^, (46) 
p=o 

where G^^\0) represents G'*-^''(A)|a=o = d^G /dX^\x=o, i.e., the derivative of order P of 
the the function G{X) at A = 0. From Eqs. (H5!) and (l46l) . the coefficient of A^ in the 
expansion of G{X) is A(P) = tr B(p)e^ ■ Therefore, we have 

A(P) ^ (47, 



It is easy to directly evaluate Eq. ( 1451) to obtain 

k 

f2 /q°° ln(l— Ae~^)a;-^cia; 



= e 

Using the expansion of 

ln(l - = - le -le--- (forlel < 1), (49) 

we obtain 

G{X) = e'''^^-^^'\ (50) 
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It can be shown that 
G(°)(0) 

G(2)(0) 



21]G'(°)(0) 
2VL 



G(-)(0) = 21](P-l)!g— 

j=o ■ ^ ' 



'(0). 



(51) 
(52) 

(53) 
(54) 



Using these recursion relations and Eqs. (HOl) and (H71) . we can then evaluate the 
expectation value of S^. In the same way, the reduced density matrix can be written as 



Ps{t) 



Pll(t) Pl2(t) 
P2l(t) P22(t) 



(55) 



where 



^ oo oo 

pii(t) = ^Y.H (I'^r^^* + H^B* + 5*7^*5 + ItI'^S*) A(Pi)A(^^) 

Pl=0 P2=0 

oo oo 

Pi2(t) = ^ E E (1^1'^^* + '^7*^^ + ^*7i?i?* + \l\^AB) AiP^)A{PS7) 



Pl=0P2=0 
P2l(i) = Puit), 



{5t 



^ oo oo 

P22(t) = ^Y.Y1 il^l^BB* + 57MP + 6*-fA*B* + \j\^AA*) A(Pi)A(^^.) 



Pl=0P2=0 



3. Results and Discussion 



3.1. Spin inversion and Rabi oscillation decay 

In this section, we present and discuss the results we obtain. We first study the dynamics 
of the central spin inversion. Figures [T](a) and[T](b) show the time evolution of (Sq) for 
different values of fl in the case of resonance (i.e., detuning e = 0). Figure [2]^a) shows 
the time evolution of (5*0) also in resonance (i.e., detuning e = 0) but with a different 
value of the system-environment coupling strength Jq from that of Fig. [I](a). For the 
parameters chosen in Figs. W^a), [11(b) and[2t^a), the driving strength g is much larger 
than the coupling strength, i.e., g ^ Jq. As a result, the self-Hamiltonian is dominant 
over the interaction with the environment. The eigenstates of the self-Hamiltonian are 
1+) = (|1) + \0))/\/2 and |— ) = (|1) — |0))/-\/2, separated by a large Rabi frequency 
g. The main influence of the environment on the central spin is to destroy the initial 
phase relation between the states |-|-) and |— ). This leads to the decay of (Sq), i.e., Rabi 
oscillation decay. From Figs. HJ^a) and[2]^a), we can see that as expected, increasing the 
value of Jo results in the increase of the decay rate of the Rabi oscillations. Apart from 
the coupling constant Jq, the important factor fl, Eq. (J26l) . also reflects the influence of 
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gt gt 

Figure 1. Time evolution of (S'q) for different values of {a) Vt — 2 and (b) VI — 30. 
The initial qubit state is IV'(O)) — |1) and other parameters are e = and Jq — Q.Qhg. 




20 40 60 80 20 40 60 80 



gt gt 

Figure 2. (a) Time evolution of (5q) for different values of 17 = 2, Jq = O.I5 (green 
sohd curve), = 2 x 10, Jq = 0.1.g/-\/lO (blue dashed curve) and = 2 x 100, 
Jo = O.lg/VlOO (red dotted curve). The initial qubit state is |'0(O)) — |1) and the 
other parameter is e = 0. (b) Time evolution of (S'q) for different values of 17 = 1, 
Jo = 0.05g (green solid curve), = 1 x 10, Jq = 0.055/\/l0 (blue dashed curve) and 
= 1 X 100, Jo = 0.055/yiOO (red dotted curve). The initial qubit state is |V'(0)) = |1) 
and the other parameter is £ = 0.5g. The inset shows the time evolution of the case of 
VL = 1 and Jo = 0.05(? in a larger time interval. 



the environment on the central spin as shown in Fig. [Hb). We can see from Figs. Hl^a), 
[Il^b) and |2t^a) that increasing the factor Vt and increasing the couphng constant Jq have 
similar effects. One can observe that the larger the value of VL is, the stronger the decay 
of the amplitude of the Rabi oscillations will be. As in the case of the spin-boson model 
discussed in Ref. [56], the central spin inversion does not oscillate around the value of 
(S'q) = 0. Its oscillations are, however, biased (shifted) a little bit toward the positive 
value of the initial (5*0(0)). With the increase of the value of Vt or Jq; this effect is 
enhanced. 

It was shown in Ref. [17] that the form and the rate of Rabi oscillation decay are 
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useful in experimentally determining the intrabath coupling strength for a broad class 
of solid-state systems [17]. This is also the case in our problem. The central spin is 
a two-level system. However, due to its interaction with the environment, the time 
evolution of the central spin inversion consists of different frequencies involved in the 
time series of Eqs. (H0ll -(H2l) as well as Eqs. (I33l)-(l35ll. The frequencies from Eq. (l35il 
can be written as 

= \l[UP2-Pi)-ef + 9'' (60) 

for a pair of integers Pi and P2. The probability distribution of the frequencies is then 

^ 00 

^(«^) = ^E^(^ + ")^(^)' (61) 
p=o 

where n = \P2 — Pi\- It is possible that other pairs of integers P[ and P2 may exist 
which correspond to the same frequency k but with m = \P2 — P[\ ^ n. In such case, 

we should add an additional probability 

^ 00 

a(«:) = -^A(P + m)A(P). (62) 
p=o 

For example, the frequencies of Eq. fl60|) can be rewritten as 

K = ^JM{P2-Pl)-{e/J^)f + g\ (63) 

If n = IP2 — Pi| is chosen to be 3, then other pairs of integers P[ and P2 with 
m = IP2 — P[\ = 17 correspond to the same frequency k for the parameters e = O.bg 
and Jo = 0.05(7 used in Fig. Mjo). Thus, summing the probability distributions with all 
possible combinations of Pi and P2 that correspond to the same frequency k leads to 
the final frequency probability distribution. 

Figure M^a) shows the frequency probability distribution for the case of e = 0, 
Jq = 0.05(7, and Q = 5. It is obvious that the main frequencies are located near g 
which is approximately the Rabi frequency. This is also the case for Figs. [I](a), [T](b) 
and [2](a) where the detuning 6 = 0. The contribution of many different frequencies, a 
consequence of interacting with the AF bath, results in the Rabi oscillation decay of 
the central spin. For the case of non-zero detuning e = 0.5g, Fig. [2](b), with Jq = 0.05g 
and = 1, shows a decay of oscillations with an approximated Rabi frequency of 
\/e^~+~g^ ~ 1-12(7. The oscillation residual amplitude approaches almost zero for a time 
period of about 65 < gt < 80 [see also the inset of Fig. [2]|^b)]. In such case, different 
frequencies interfere with each other and produce the zero amplitude variation period, a 
behavior that is called collapse in quantum optics [58]. At longer times, the oscillation 
amplitude revives. This is shown clearly in the inset of Fig. Mjo). Reference [17] reported 
an experimental observation of a similar collapse and revival phenomenon for driven spin 
oscillations, though the spin bath and the corresponding interactions are different from 
those discussed here. 

The environmental conditions affect the dynamics of the central spin and the 
completeness of the collapse. Figure Et^b) shows the frequency probability distribution 
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K/g K/g 

Figure 3. (a) Probability distribution of frequencies for e = 0, = 5, and Jq = Q.Qbg. 
The inset shows the comparison between the probability distribution of frequency 
and that of the case of increasing N 10 times (i.e., = 5 x 10, Jq = O.OSg/VTO) 
with other parameters unchanged, (b) Probability distribution of frequencies for the 
parameters used in Fig.[2jb). The inset shows the comparison between the probability 
distribution of frequencies and that of the case of increasing N 10 times (i.e., $7 = 1 x 10, 
Jo = O.OSg/vTO) with other parameters unchanged. 

for the central spin inversion with parameters used in Fig. Et^b). The distribution has 
a left- hand- side cutoff at Umin = 9 and the center of the distribution is located at 
Kc = a/^^ + g^- One can easily understand this from the frequency relation Eq. ( l60l) . 
It is also obvious from Eq. (!60|) that the center of the distribution shifts toward larger 
frequencies with the increase of the detuning [this can also be seen from the comparison 
between Figs. [21(a) and [31(b)]. If the detuning e exists, it is possible that the shape of the 
frequency probability distribution becomes approximately Gaussian, which then results 
in well-defined collapse and revival behavior regions [SB]. If the detuning is zero, only 
a half side of the distribution exists [see Fig. [3]^a)] and the central spin inversion will 
never show the (complete) collapse and revival behaviors. The frequency probability 
distribution is determined also by the coupling strength Jq and the important factor 
Q. With the increase of Jq or Q, the width of frequency distribution increases and the 
probability decreases. As a result, the decay of the Rabi oscillations is enhanced. 

We note that the results presented here depend, of course, on the number of 
environment atoms N in each sublattice. The influence of the number of atoms N 
in each sublattice on the dynamics of the central spin comes from two respects. One 
is the scaled coupling constant Jq = Jq//\/N which is proportional to 1/a//V. The 
other is the factor Q, Eq. ( l26l) . which is proportional to A^. We plot in Fig. [2]^a) 
and Fig. Mjo) the time evolutions of (Sq) in (blue) dashed curve and (red) dotted 
curve corresponding, respectively, to increasing by 10 and 100 times but keeping 
other parameters unchanged. The results show that the two curves in dashed and 
dotted lines almost coincide with each other. If we increase A^ by 200 times, then the 
resultant evolution and that of increasing A^ by 100 times become indistinguishable and 
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Figure 4. (a) Entropy for different initial qubit states: ["(/"(O)) = (|v?i) + |v52))/\/2 
(pink dot-dashed curve), |V'(0)) = |1) (blue dashed curve), |'(/'(0)) = |<(5i) or |V-'(0)) = 
|</72) (red solid curve). Other parameters are e = g, Jq ^ O.Olg, ~ 20. (b) Entropy 
E{t) for (i) pink dot-dashed curve: Jo = 0.03g, e = 0, (ii) green dashed curve: 
Jo — 0.01.9, e 0, (iii) red sohd curve: Jq = O.Olg, e = g, and (iv) blue dotted 
curve: Jq = 0.01(7, e — 3g. The qubit initial state is |'0(O)) ~ \(pi) or \ip{0)) = \^2), 
and other parameter is f2 = 20. 



converge toward the same evolution. In other words, a nontrivial finite hmit exists in the 
thermodynamic hmit of — ?■ oo. This has been shown analytically in Refs. [261 123 [2H] 
for a spin coupled to an AF environment without an external ac driving field. We show 
numerically here that this N oo limit also exists in the driven spin case. One can also 
notice that the difference between the original curve (in solid line) and the two curves 
(in dashed and dotted lines, respectively) in Fig. [2] is very small in the all time interval 
shown. This difference in Fig. EJ^b) with e ^ is also small but is slightly larger at 
large times. These indicate that the results presented here are very close to the results 
in the thermodynamic limit. If we also increase the AF atom number A^ in Jq and Q 
simultaneously but keep other parameters unchanged, then we can find that the number 
of the contributed frequencies in the probability distribution of frequencies becomes 
larger (i.e., the separation between two adjacent distributed frequencies becomes smaller 
or denser) and the value of the distribution becomes smaller [see the insets in Figs. M^a) 
and|3]^b)]. But the width and the shape of the frequency probability distribution remain 
almost the same as before [see the insets in Figs. EJ^a) and Et^b)]. The fact that the 
frequency distribution does not change much results in the very similar dynamics of the 
driven spin (see Fig. |2]) when we increase A^. The existence of nontrivial finite limits as 
N ^ oo and the small difference between the original dynamical evolutions with their 
counterparts in the thermodynamic limit are true for all the dynamic variables {{Sq) 
and von Neumann entropy) presented in this article. 
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3.2. von Neumann entropy 




(64) 



After diagonalizing the reduced density matrix ps{t), we obtain the von Neumann 



fc=i 

where Pkit) are the eigenvalues of the reduced density matrix ps(t) and can be expressed 
as 



We investigate the von Neumann entropy of the qubit system in Fig. |H In Fig. Hl^a), for 
the initial state IV'(O)) = the entropy E{t) oscillatorily approaches its maximal value 
of 0.4228 (in blue dashed curve). If the initial state is one of the eigenstates of the qubit 
self-Hamiltonian, i.e., \ipi) or \(f2), the von Neumann entropy E{t) remains small with 
time and approaches a very small value in the long-time limit (in red solid line). In the 
opposite, if we choose the superposition state |'?/'(0)) = + \ip2))/V^ as the initial 
state, results show that the von Neumann entropy E{t) increases monotonously to the 
maximal value of In 2 = 0.6971 (in pink dot-dashed line), representing a completely 
mixed state. That is to say, this initial state gets maximally entangled with the 
environment in the course of the evolution. This initial superposition state is thus 
frangible to the influence of the environment and should be avoided being used in the 
quantum information processing in the presence of the AF environment. In all the cases, 
the entropy E(t) increases with the increase of the factor fl and the coupling constant 
strength Jq. 

We discuss next the pointer states of the qubit from the perspective of the von 
Neumann entropy [591 EOl EI]- For 9 = 0, i-e., no driving field, it is obvious that the 
eigenstates of the self-Hamiltonian, ^qSq, are the pointer states as the self-Hamiltonian 
commutes with the interaction Hamiltonian with the environment. In this case, a 
generic quantum state decays, after the decoherence time, into a mixture of pointer 
states. Thus the decoherence behavior can be described effectively by the off-diagonal 
elements (decoherence factor) of the reduced density matrix in the pointer state basis 
of the system. This is exactly the cases studied in Refs. [261 EZ], where the off-diagonal 
elements (decoherence factor) vanish in the long-time limit in the eigenstate basis of the 
self-Hamiltonian or of the interaction Hamiltonian with the environment regardless of 
the different initial states. If the driving field strength 7^ 0, it is difficult to identify the 
exact pointer states of the system (especially in the non-perturbative and non-Markovian 
regime) as the self-Hamiltonian does not commute with the interaction Hamiltonian with 
the bath. Pointer states may be defined as the states which become minimally entangled 
with the environment in the course of their evolution [59l|60]. An operational definition 



entropy 



2 




(65) 



Pnjt) + P22{t) ± VlPujt) - p22{tW + 4pi2(t)p2l(t) 

2 



(66) 
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in terms of the von Neumann entropy, introduced in Refs. [5^1 EOl EI], is that pointer 
states are obtained by minimizing the von Neumann entropy over the initial state |^/'(0)) 
and requiring that the answer be robust when varying the time t. In general situations, 
pointer states result from the interplay between self-evolution and interaction with the 
environment, and thus their dynamical selection by the environment are complicated. 
For g » Jq or/and e » Jq, the pointer states turn out (approximately) to be the 
eigenstates of the self-Hamiltonian [211 ES]- Here the eigenstates \ipi) and \ip2) of the 
self-Hamiltonian Hs = bSq + gS^ defined in Eq. ([5]) can be written as 

(67) 
(68) 





where 



Uii = ^ . (69) 



e'^ + g'^ — eJe^ + g^ 



V2 
1 

7! 



2 



^2 _^ g2 _ S^S2 _j_ g 

U21 = -4 , ^ (71) 

+ (^2 _|_ £^52 _|_ g2 



" " 7^/ 2 2 r-2 -2 ^ ' 

Thus for the Hamiltonian parameters chosen in Fig. Hl^a) which is in a regime where 
the self-Hamiltonian of the system dominates, it is expected that the pointer states 
are very close to the eigenstates Iv^i) and 1^92) [Ml ISS]- By changing different initial 
states near the eigenstates \^\) or \^f2) with other parameters fixed, it is found that the 
eigenstate or |(^2) produces the minimum entropy increase with time (the red solid 
line). This suggests that the eigenstates |</9i) and 1^92) are dynamically selected by the 
environment as the preferred pointer states in the self-Hamiltonian dominant regime. 
Moreover, when the energy parameters g and e of the self-Hamiltonian of the central 
spin system are increased further with respect to the system-bath coupling strength Jq 
(i.e., the self-Hamiltonian dominates even more), the values of entropy evolution are 
lower (the green dashed curve is lower than the pink dot-dashed curve in the zero- 
detuning case, and the blue dotted curve that is close to zero is lower than the red solid 
curve in the detuning case) as shown in Fig. 111(b). We have also checked that the four 
curves with initial state or |(/?2) in Fig. Hl^b) are, respectively, the minimum entropy 
curves in the course of evolution among those entropy evolution curves with initial 
states varied near the eigenstate of \ip\) or but with other parameters fixed. This 
confirms furthermore again that in the regime where the self-Hamiltonian of the system 
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dominates (with or without detuning), the eigenstates of the system self-Hamiltonian 
emerge as the preferred pointer states. 

4. Conclusions 

We investigate the decoherence of a qubit in an AF environment, in the presence of a 
driving field. The difficulty of our problem hes in the fact that the self-Hamiltonian does 
not commute with the interaction Hamiltonian and the internal dynamics (coupling) 
of the spin bath are taken into consideration at the same time. The spin-wave 
approximation is used to map the spin operators of the AF environment onto bosonic 
operators in the low-temperature and low-excitation limit. Then the resultant model 
is solved exactly, even in the case of mult i- environment modes and finite environment 
temperatures. Our approach includes the environment dynamics, the qubit dynamics, 
and the quantum correlations between them. The influence of the AF environment 
on the qubit is found to depend on two factors, i.e., the coupling constant Jq and a 
dimensionless factor Q = ^^^'^3/2 js ■ Increasing the two factors, the decay time of the 
Rabi oscillations becomes shorter, the decoherence of the qubit is enhanced, and the 
qubit state becomes more mixed with time. The time evolution of the Rabi oscillations 
also depends on the detuning between the driving frequency and the qubit Larmor 
frequency. If the detuning exists, the Rabi oscillations may show a behavior of collapses 
and revivals; however, if the detuning is zero, such a behavior will not appears. This can 
be understood in terms of the weighted frequency distribution investigated here. Also 
the decoherence and the pointer states of the qubit are discussed from the perspective of 
the von Neumann entropy. It is found that the eigenstates of the qubit self-Hamiltonian 
are dynamically selected by the environment as the preferred pointer states in the weak 
system-environment coupling limit (or in the self-Hamiltonian dominant regime). 

Acknowledgments 

H.S.G. would hke to acknowledge support from the National Science Council, Taiwan, 
under Grant No. 97-21 12-M-002-012-MY3, support from the Frontier and Innovative 
Research Program of the National Taiwan University under Grants No. 99R80869 and 

No. 99R80871, and support from the focus group program of the National Center for 
Theoretical Sciences, Taiwan. H.S.G. is also grateful to the National Center for High- 
performance Computing, Taiwan, for computer time and facilities. X.Z.Y. acknowledges 
support from the National Science Foundation of China under Grant No. 10874117. 
X.Z.Y. also acknowledges the support by the National Science Foundation (Grant No. 
11091240282) to do research in ICTP (SMR2154), and thanks ICTP and APCTP for 
their invitations as a visiting professor and for their hospitality. 

References 

[1] B. E. Kane, Nature (London) 393, 133 (1998). 



Dynamics of a driven spin coupled to an antiferromagnetic spin hath 



18 



[2] H.-S. Goan, Int. J. Quantum Inf. 3, 27 (2005); L. C. L. HoUcnbcrg ct al., Phys. Rev. B 74, 045311 

(2006) ; C. J. Wellard et al, ibid. 68, 195209 (2003); B. Koiller et al, Phys. Rev. Lett. 88, 
027903 (2001); L. M. Kettle et al., Phys. Rev. B 73, 115205 (2006). 

[3] C. D. Hill and H.-S. Goan, Phys. Rev. A 68, 012321 (2003). G. D. Hill and H.-S. Goan, Phys. 

Rev. A 70, 022310 (2004); L. M. Kettle et al., Phys. Rev. B 68, 075317 (2003). 
[4] D.-B. Tsai, P.-W. Chen, and H.-S. Goan, Phys. Rev. A 79, 060306(R) (2009). C. D. Hill et al., 

Phys. Rev. B 72, 045350 (2005). 
[5] G. Burkard, D. Loss, and D. P. DiVincenzo, Phys. Rev. B 59, 2070 (1999). 
[6] R. Rausscndorf and H. J. Briegel, Phys. Rev. Lett. 86, 5188 (2001). 
[7] X. Hu and S. Das Sarma, Phys. Rev. Lett. 96, 100501 (2006). 
[8] D. Loss and D. P. DiVincenzo, Phys. Rev. A 57, 120 (1998). 

[9] V. Cerletti, W. A. Coish, O. Gywat, and D. Loss, Nanotechnology 16 R27 (2005). 
[10] M. Borhani, V. N. Golovach, and D. Loss, Phys. Rev. B 73, 155311 (2006). 
[11] A. Tackeuchi, T. Kuroda, K. Yamaguchi, Y. Nakata, N. Yokoyama, and T. Takagahara, Physica 
E 32, 423 (2006). 

[12] L. Gywinski, W. M. Witzcl, and S. Das Sarma, Phys. Rev. B 79, 245314 (2009). 

[13] W. Zhang, N. P. Konstantinidis, V. V. Dobrovitski, B. N. Harmon, L. F. Santos, and L. Viola, 

Phys. Rev. B 77, 125336 (2008). 
[14] W. Zhang, V. V. Dobrovitski, K. A. Al-Hassanieh, E. Dagotto, and B. N. Harmon, Phys. Rev. B 

74, 205313 (2006). 

[15] W. Zhang, V. V. Dobrovitski, L. F. Santos, L. Viola, and B. N. Harmon, Phys. Rev. B 75, 201302 

(2007) . 

[16] J. Lages, V. V. Dobrovitski, M. I. Katsnelson, H. A. De Raedt, and B. N. Harmon, Phys. Rev. E 
72, 026225 (2005). 

[17] V. V. Dobrovitski, A. E. Feiguin, R. Hanson, and D. D. Awschalom, Phys. Rev. Lett. 102, 237601 

(2009). 

[18] S. Takahashi, R. Hanson, J. van Tol, M. S. Sherwin, and D. D. Awschalom, Phys. Rev. Lett. 

101, 047601 (2008). 
[19] W. Yang and R.-B. Liu, Phys. Rev. B 78, 085315 (2008). 
[20] W. Yang and R.-B. Liu, Phys. Rev. B 79, 115320 (2009). 
[21] W. M. Witzel and S. Das Sarma, Phys. Rev. B 74, 035322 (2006). 
[22] V. V. Dobrovitski and H. A. De Raedt, Phys. Rev. E 67, 056702 (2003). 

[23] K. A. Al-Hassanieh, V. V. Dobrovitski, E. Dagotto, and B. N. Harmon, Phys. Rev. Lett. 97, 

037204 (2006). 

[24] F. M. Cucchictti, J. P. Paz, and W. H. Zurek, Phys. Rev. A 72, 052113 (2005). 

[25] J. Dziarmaga, Phys. Rev. B 71, 054516 (2005). 

[26] X. Z. Yuan and K. D. Zhu, Europhys. Lett. 69, 868 (2005). 

[27] X. Z. Yuan, H. S. Goan, and K. D. Zhu, New J. Phys. 9, 219 (2007). 

[28] X. Z. Yuan, H.-S. Goan, and K. D. Zhu, Phys. Rev. A 81, 034102 (2010). 

[29] E. Novais, A. H. Castro Ncto, L. Borda, I. Affleck, G. Zarand, Phys. Rev. B 72, 014417 (2005). 

[30] H. P. Brcucr, Phys. Rev. A 69, 022115 (2004). 

[31] A. Button and S. Bose, Phys. Rev. A 69, 042312 (2004). 

[32] M. Lucamarini, S. PaganeUi, S. Mancini, Phys. Rev. A 69, 062308 (2004). 

[33] H. P. Brcucr, D. Burgarth, and F. Petruccionc, Phys. Rev. B 70, 045323 (2004). 

[34] Y. Hamdouni, M. Fannes, and F. Petruccionc, Phys. Rev. B 73, 245323 (2006). 

[35] X. Z. Yuan, H.-S. Goan, K. D. Zhu, Phys. Rev. B 75, (2007). 

[36] H. T. Quan, Z. Song, X. F. Liu, P. Zanardi, C. P. Sun, Phys. Rev. Lett. 96, 140604 (2006). 
[37] F. M. Cucchictti, S. F. Vidal, J. P. Paz, Phys. Rev. A 75, 032337 (2007). 

[38] D. Rossini, T. Calarco, V. Giovannetti, S. Montangero, R. Fazio,, Phys. Rev. A 75, 032333 (2007). 

[39] H. Krovi, O. Oreshkov, M. Ryazanov, and D. A. Lidar, Phys. Rev. A 76, 052117 (2007). 

[40] E. Ferraro, H.-P. Breuer, A. Napolil, M. A. Jivulescu, and A. Messina Phys. Rev. B 78, 064309 



Dynamics of a driven spin coupled to an antiferromagnetic spin hath 



19 



(2008) 

[41] N. Arshed, A. H. Toor, and D. A. Lidar, Phys. Rev. A 81, 062353 (2010) 

[42] J. R. Petta, A. C. Johnson, J. M. Taylor, E. A. Laird, A. Yacoby, M. D. Lukin, C. M. Marcus, 

M. P. Hanson, and A. C. Gossard , Science 309, 2180 (2005). 
[43] M. H. Mikkelsen, J. Berezovsky, N. G. Stoltz, L. A. Coldren, and D. D. Awschalom, Nature Phys. 

3, 770 (2007). 

[44] F. H. L. Koppens, C. Buizert, K. J. Tielrooij, I. T. Vink, K. C. Nowack, T. Meunier, L. P. 

Kouwenhoven, and L. M. K. Vandersypen, Nature (London) 442, 766 (2006). 
[45] B. C. Stipe, H. J. Mamin, C. S. Yannoni, T. D. Stowe, T. W. Kenny, and D. Rugar, Phys. Rev. 

Lett. 87, 277602 (2001). 

[46] F. H. L. Koppens, D. Klauser, W. A. Coish, K. C. Nowack, L. P. Kouwenhoven, D. Loss, and L. 

M. K. Vandersypen, Phys. Rev. Lett. 99, 106803 (2007). 
[47] R. Hanson, V. V. Dobrovitski, A. E. Feiguin, O. Gywat, D. D. Awschalom, Science 320, 352 

(2008). 

[48] D. D. Bhaktavatsala Rao, Phys. Rev. A 76, 042312 (2007). 

[49] M. Lucamarini, S. PaganclU, and S. Mancini, Phys. Rev. A 69, 062308 (2004). 

[50] A. Lupa§cu, E. F. C. Driessen, L. Roschier, C. J. P. M. Harmans, and J. E. Mooij, Phys. Rev. 

Lett. 96, 127003 (2006). 
[51] C. Kittcl, Quantum Theory of Solids, (John Wiley and Sons, Inc., New York, 1963). 
[52] H. P. Breuer, D. Burgarth, and F. Petruccione, Phys. Rev. B 70, 045323 (2004). 
[53] M. Prasca, Ann. Phys. N.Y. 313, 26 (2004). 

[54] K. Yosida, Theory of Magnetism, (Springer series in solid-state sciences, V. 122). 
[55] A. J. Lcggctt, S. Chakravarty, A. T. Dorsey, M. P. A. Fisher, A. Garg and W. Zwerger, Rev. 
Mod. Phys. 59, 1 (1987) 

[56] E. K. Irish, J. Gea-Banacloche, I. Martin, and K. C. Schwab, Phys. Rev. B 72, 195410 (2005). 

[57] H. Xiong, S. Liu, G. Huang, and Z. Xu, Phys. Rev. A 65, 033609 (2002). 

[58] M. O. Scully, M. S. Zubairy, Quantum Optics, (Cambridge University Press, Cambridge, 1997). 

[59] J. P. Paz, and W. H. Zurek, Phys. Rev. Lett. 82, 5181 (1999). 

[60] W. H. Zurek, S. Habib, and J. P. Paz, Phys. Rev. Lett. 70, 1187 (1993). 

[61] W. H. Zurek, Rev. Mod. Phys. 75, 715 (2003). 



